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We show in the case of a superconducting Nb ladder that a mesoscopic superconductor typ- 
ically exhibits magnetoresistance oscillations whose amplitude and temperature dependence are 
different from those stemming from the Little-Parks effect. We demonstrate that these large resis- 
tance oscillations (as well as the monotonic background on which they are superimposed) are due 
to current- excited moving vortices, where the applied current in competition with the oscillating 
Meissner currents imposes/removes the barriers for vortex motion in increasing magnetic field. Due 
to the ever present current in transport measurements, this effect should be considered in parallel 
with the Little-Parks effect in low-T c samples, as well as with recently proposed thermal activation 
of dissipative vortex-antivortex pairs in high-T c samples. 

PACS numbers: 74.78.Na, 74.25. F-, 74.25.N-, 74.40.Gh 



The circulation of electrons and Cooper pairs in dou- 
bly and multiply connected normal and superconducting 
materials in the presence of a magnetic field produces re- 
spectively the Aharonov-Bohm [T] and the Little-Parks 
(LP) effect 0. These effects are sensitive probes for the 
proposed incoherent Cooper pairing in the pseudogap [5] 
and for insulating phases [I] competing with supercon- 
ductivity in copper oxide and conventional superconduc- 
tors, for the relative contribution of individual bands in 
two-band superconductors [5] and for the interplay of su- 
perconductivity and magnetism in hybrid structures [BJ. 
They are expected to shed light on the d-wave symmetry 
[7] and the strip structures [8] in high temperature su- 
perconductors. Magnetoresistance oscillations have been 
widely observed in doubly and multiply connected meso- 
scopic superconductors [9TfTB] and have mostly been at- 
tributed to the LP effect, which is related to the periodic 
suppression of the critical temperature T c due to fluxoid 
quantization (TTl [T8| . However, analysis of magnetore- 
sistance oscillations in recent experiments on high-T c su- 
perconducting loops (in a specifically designed network) , 
showed that neither the amplitude of the oscillations, nor 
their temperature dependence could be explained by the 
LP effect [15] . In order to describe the experimental data, 
the authors used the fluxoid dynamics model [16] , accord- 
ing to which thermally excited vortices/antivortices move 
and interact with the magnetic field-induced persistent 
currents. 

Here, we present results of numerical simulations and 
transport measurements, which reveal a new origin for 
the magnetoresistance oscillations in mesoscopic super- 
conductors, namely, the motion of current-excited vor- 
tices. These vortices are not excited by thermal fluctu- 
ations in contrast to the case of Refs. [HI [16]. Instead, 
the applied current interacts with existing Meissner cur- 
rents in the sample and tunes the barrier for vortex en- 



try and exit while simultaneously driving the vortices. 
In addition, the ubiquitous monotonic background of the 
magnetoresistance curves is also found to be due to mov- 
ing vortices. Since an external current is inherent to all 
transport measurements, the proposed mechanism needs 
to be considered in the studies of singly (stripes), dou- 
bly (loops) or multiply (ladders and wire networks) con- 
nected mesoscopic superconductors. 

In what follows, we demonstrate our findings for a Nb 
ladder structure (i.e., a thin strip with a chain of holes, 
see Fig. [T]), which is a transition geometry from a fi- 
nite ring structures to an infinite network [19) . In this 
geometry, the electrical contact issues on the resistance, 
known to occur in individual loops (20) . are suppressed. 
The choice of material and geometry was made to maxi- 
mize the larger characteristic length scales, to avoid the 
formation of "effective" weak links, which can result in 
direct superconducting-insulator transitions |21j and con- 
sequently large magnetoresistance oscillations in a broad 
range of temperatures and magnetic fields |22j . Finally, 
thermal fluctuations in Nb are far less important than in 
high-T c samples, which is also relevant for the analysis of 
our findings in comparison with Ref . [15] . 

Theory. - For the theoretical study of the fluxoid 




FIG. 1: (color online) An oblique view of the system: a su- 
perconducting strip (of width w = a + 25, thickness f < A,^ 
and periodically long in the y-direction) with a chain of holes 
(size a and period d = a + 5) under applied dc current I and 
a perpendicular magnetic field H. 
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quantization and dynamics, we used the generalized time- 
dependent Ginzburg-Landau (GL) equations |23j : 
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with units: the coherence length £(0) for the distance, 
iGi(0)=7rfi,/8fcsT'cM for time, T c for temperature, and 
<^o/27r£(0) for the vector potential A. Order param- 
eter ip is scaled to its value for zero magnetic field 
and temperature, and the current density is in units of 
jo = c n ft/2e£GL(0)£(0)- The material parameters 7 = 20 
and u = 5.79 follow from microscopic theory [23]. x(r,t) 
is the random fluctuating term [24] . To account for heat- 
ing effects, we couple Eqs. (1-2) to the heat transfer 
equation (see Ref. [20] for details): 
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where To is the bath temperature. Here we use ^=0.03, 
C=0.06 and rj = 2-10~ 4 , corresponding to an intermediate 
heat removal to the substrate [5D]. The above equations 
are solved sclf-consistently using the semi-implicit Crank- 
Nicholson algorithm [25) with periodic boundary condi- 
tions in the y-direction and Neumann boundary condi- 
tions at all sample edges. The transport current is in- 
troduced via the boundary condition for the vector po- 
tential, rotA\ z (x—0,w)=HzLHi, where H]=2wl /c is the 
magnetic field induced by the current /. 

As a representative example, we studied the response 
to an external magnetic field of a superconducting strip 
of width w = 40£(0) and with holes of size a — 20£(0). 
Figure [2|a) shows the time averaged energy F (see Ref. 
[2"6] for the definition of F) of this system as a function 
of magnetic field for different values of the applied dc 
current I. For small current values (curve 1), the vortex 
entries are shifted to lower magnetic field compared to the 
case of no applied current (not shown here) . Still, vortex 
entry remains a first-order transition, i.e., one observes a 
sharp drop in the energy at the transition field (indicated 
by a solid black arrow in Fig. [2]ja)). However, for the 
same current at larger magnetic fields we surprisingly 
found that vortex entry is characterized by a continuous 
change in energy, manifested by a second-order transition 
(dashed blue arrow). It is exactly these transitions that 
give rise to the large voltage signal V — — J dA/dt-dl, as 
shown in Fig. [2|b) (see curve 1). Such vortex entries are 
found at lower fields for larger applied current (curves 2 in 
Fig. |2ja)). With further increasing /, the magnetic fields 
at which voltage oscillations are located do not change, 




FIG. 2: (color online) Time-averaged free energy F (a) (in 
units of Fo = HcVo/8ir with Vb the unit cell volume) and 
voltage V (b) (in units of normal state voltage V n ) as a func- 
tion of external flux (f> (calculated over the unit cell area 
S — w x (a + 8)) for different values of applied current 7. 
The width of the simulated ladder is w — 40£(0) and the size 
of the holes is a = 20£(0). The temperature is T = 0.98T C 
and the GL parameter is « = 10. 



but oscillations become broader (curves 3 and 4 in Fig. 
[2jb)). Note that oscillations in our curves are not exactly 
periodic with field and larger flux is needed to observe 
the first voltage peak (a typical effect for finite-walled 
ring structures, see e.g. Ref. [TT]). 

To gain further insight into the process leading to mag- 
netoresistance peaks, we plot the temporal voltage signal 
for I/{j w) = 0.04 and <j> = 0.6400 in Fig. [3] At this field 
(i.e., just after the maximum of the first voltage peak), 
the output voltage oscillates periodically in time with a 
global minimum (point 1 in Fig. |3| corresponding to one 
vortex trapped inside the hole (inset 1). With time, this 
vortex exits the sample (inset 2) exhibiting a maximum 
in the V(t) characteristics (point 2), and the system is 
free of vortices (inset 3). Immediately after, a new vor- 
tex penetrates the sample (inset 4) and the entire vortex 
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FIG. 3: (color online) Voltage vs. time response for the sam- 
ple of Fig. [5] for cj) = O.640o and I/(j w) = 0.04. Insets show 
snapshots of \ip\ 2 at time intervals indicated by open circles 
on the V(t) curves. 
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FIG. 4: (color online) (a,b) Time evolution of the free en- 
ergy for the ladder considered in Fig. [5] at <j> = 0.54<^o (a) 
and (j> — 0.71(?!)o (b). Arrows indicate the penetration AF P 
and expulsion AF e barriers, (c) Time-averaged voltage (solid 
curve), AF P (filled circles) and AF e (open circles) as a func- 
tion of external flux. The applied current is I/(jow) = 0.04. 



entry-exit sequence repeats (see Ref. [27] for an animated 
data for this process). 

We stress once more that vortices in our study cannot 
be induced by thermal fluctuations, making our predic- 
tions different from ones in Refs. [HI [16]. Note also that 
no voltage/resistance oscillations with time are expected 
from the LP effect, which is also different from our find- 
ings. The origin of vortex excitation in our system is the 
interplay of currents. Namely, applied magnetic field in- 
duces reactive screening (and circulating) currents in the 
superconductor, j s . Therefore, applied dc current always 
enhances the supercurrent j s on one side of the sample, 
while suppressing j s on the other side. This directly af- 
fects the barriers for vortex entry and exit, which are 
indicated in Figs. [4^a,b) , where we plotted the time evo- 
lution of the free energy of the system. When the total 
current j t = j s + j reaches to its critical value [28] in 
one arm of the ladder, a vortex nucleates at the edge of 
the sample in spite of the finite energy barrier AF p (see 
filled circles in Fig. |4|c)). At that field the barrier for 
the vortex expulsion AF e on the opposite side of the sam- 
ple is strongly suppressed (see open circles in Fig. (4^c)), 
since jt — j s — j- Therefore, the vortex is driven through 
the system resulting in periodic-in-time voltage signal as 
shown in Fig. [3] With further increasing the applied 
field, the barrier for vortex exit increases (open circles in 
Fig. |4]^c)). Consequently, the vortex remains trapped in 
the sample and voltage becomes zero. Such modulation 
of the barrier for vortex entry and exit is observed for the 
transition between different vortex states (see the second 
voltage peak in Fig. |4^c)). 

Thus, we conclude that each voltage (and magnetore- 



sistance) peaks as a function of magnetic field corre- 
sponds to the moving vortices in the system, the motion 
of which is manipulated by the interplay of applied and 
Meissner currents. Note however that at higher driving 
current or for larger vorticity, it is no longer possible 
to stabilize the stationary vortex state and vortices al- 
ways remain in motion, creating a background voltage 
(characteristic of a resistive state) on which the peaks 
for every additional vortex entry are superimposed (see 
curve 4 in Fig. gb)). For this case, animated data 
for the evolution of the superconducting condensate with 
time, together with the corresponding voltage signal, are 
shown in the Supplementary Material [27]. The back- 
ground on which the magnetoresistance oscillations are 
superimposed has been puzzling scientists since the orig- 
inal experiment of Little and Parks [2] . Several possible 
explanations have been put forth, including the Meiss- 
ner effect, misalignment of the external field, and differ- 
ent single-particle excitation energies within an electron 
pair. Recently Sochnikov et al, proposed an alternative 
mechanism which is based on thermally excited vortices 
and antivortices [15] . In our case, the background stems 
solely from the continuous, dissipative vortex motion. 

Experiment. - In order to confirm our theoretical pre- 
diction, we performed magnetoresistance measurements 
on 100 nm thick Nb strips with a chain of square holes 
patterned with focused-ion-beam milling (see Ref. [5S] 
for the details of sample fabrication and measurements). 
We considered Nb strips with width w=385 nm, hole size 
a=120 nm and period g?=300 nm, 385 nm and 800 nm. 
Magnetoresistance of the samples was obtained by four- 
probe dc measurements using a constant current mode. 
Here we present the results for the sample with inter- 
hole spacing d=385 nm (see the lower inset in Fig. [5^b)), 
which shows a broad superconducting/normal transition 
with T c ps 7.725 K (using the criterion of 90 % of the 
normal state resistance, see Fig. [5](b)) . Zero temper- 
ature coherence length £(0) and the penetration depth 
A(0) were estimated to be 10 nm and 200 nm, respec- 
tively. 

Figure [5]ja) shows experimentally measured magne- 
toresistance curves of the sample at different temper- 
atures (see Supplemental Material 30J for the results 
obtained for the sample with d=300 nm). Because of 
the small bias current (7=1 fiA and current density 
j = 2.5 ■ 10 3 A-cm~ 2 ), no resistance oscillations are ob- 
served until higher temperatures (T = 7.3 K), where a 
first resistance peak appears at around 1400 Oe. Figure 
[5jc) shows the R{H) curves obtained at T = 7.42 K for 
different values of the biased current. As delineated in 
this figure, the amplitude of the resistance oscillations 
is strongly affected by the bias current: the oscillations 
disappear both at larger and smaller currents. However, 
the position of the resistance peaks is independent of the 
applied current, as we predicted in our numerical simu- 
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FIG. 5: (color online) (a) Measured resistance (in units of 
normal state resistance R„ = 27.9 Q) vs. magnetic field for 
a Nb superconducting strip (of width w = 385 nm) with a 
chain of holes (of size a = 120 nm and period d = 385 nm), 
in applied dc current 1 = 1 /xA, and at different temperatures 
(ranging from T = 7.35 K to T = 7.65 K with step 50 mK). 
Solid curves show theoretical results for the dimensions of the 
experimental sample and for I/(jow) = 0.05. (b) R(T) curve 
of the sample at zero magnetic field for 1 = 1 /iA. Lower 
inset shows the SEM image of the sample. Upper inset shows 
the amplitude of the first resistance peak, AR, as a function 
of temperature (solid dots), compared to the amplitude of 
the resistance oscillations predicted by the Little-Parks effect, 
AR = 0.U(dR/dT)T c (£(0)/(w + o)/4) 2 [3 [15] (open dots). 
Triangles present the results of the GL simulations, (c) R(H) 
curves of the sample obtained at T = 7.42 K for different 
values of the bias current. 



lations (see Fig. |2jb)). 

We now analyze the amplitude of the resistance oscilla- 
tions as a function of temperature, shown in the top inset 
of Fig. [5^b) by solid dots for the first peak of the R(H) 
curves of Fig. [5]Ja). Open dots in this figure show the 
expected amplitude of the LP oscillations due to periodic 
changes in the superconducting transition temperature: 



AR = O.U(dR/dT)T c (Z(0)/r eff ) 2 with r eff = (w + a)/4 
[21 US] . The LP mechanism predicts resistance oscilla- 
tions at higher temperatures, T > 7.5 K, and can survive 
until very close to T c . However, we observed resistance 
peaks already starting at T = 7.35 K, and which disap- 
pear well below T c . Moreover, the amplitude of the oscil- 
lations reported here is much larger than the ones orig- 
inating from the LP effect, in spite of the fact that the 
LP effect is pronounced in low-T c superconductors due 
to their relatively large coherence length. On the other 
hand, current-driven moving vortices can account for the 
observed large magnetoresistance. The solid curves in 
Fig. [5]ja) are results of theoretical simulations for the 
parameters of the experimental sample. We found very 
good agreement with the period of the voltage oscilla- 
tions, while we observed a weak difference in the ampli- 
tude of the oscillations - likely due to the presence of 
superconducting current and voltage leads in the exper- 
iment (see triangles in the top inset of Fig. [5^b)). In 
simulations, the damage to the sample near the (hole) 
edges due to the ion bombardment process (the bright 
regions in the SEM image of the sample at the sample 
edges) was modeled by the local suppression of T c of ap- 
prox. 4 % in these regions (~ 20 nm width). We also 
used the parameter u = 1 (see Eq. (1)) to best fit the 
experimental data, which indicates a larger electric field 
penetration depth Ie in our experimental samples (since 
u ~ lg 2 ) [21] ■ The interpretation of the voltage oscilla- 



tions vs. magnetic field remains the same as in Figs. 2][4 
confirming their origin to be the moving vortices. Here 
we must emphasize that due to the small width of the 
rings (i.e., (w — a)/2 < 3£(T)) in the considered tem- 
perature range, the model of thermally created vortex- 
antivortex pairs (Ref. |15j and references therein) is not 
applicable to our experimental sample. 

In summary, we have demonstrated that the Little- 
Parks effect, a landmark manifestation of flux quantiza- 
tion, is not sufficient to fully describe the size and temper- 
ature range of typical magnetoresistance oscillations ob- 
served in small superconductors. Viable interpretations 
of such oscillations have been offered in the past by Parks 
and Mochel [32] , Anderson and Dayem [33] , Sochnikov et 
al. [T5] and Baturina et al. [55] . We however showed that 
in low-T c mesoscopic samples large magneto-resistance 
oscillations can stem solely from intermittent vortex nu- 
cleation, motion, and stabilization by the interplay of 
the driving current and the persistent Meissner currents 
in the sample. Furthermore, continuously moving (dissi- 
pative) vortices contribute to the background of the mag- 
netoresistance curves, and their contribution grows with 
increasing field or applied current - an effect not covered 
by competing explanations. We therefore suggest that, 
due to ever present current in transport measurements, 
the contribution of current-induced vortices to dissipa- 
tion in patterned superconductors can never be entirely 
disregarded. 
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